Large-scale effects on meso-scale modeling for 

scalar transport 



M. Cencini^'^, A. Mazzino^, S. Musacchio^''^ and A. Vulpiani^'^'^ 

1 ISC- CNR, Via dei Taurini 19, 1-00185 Roma, Italy 

^ INFM-SMC c/o Dipartimento di Fisica, Universitd "La Sapienza" , 
P.zzle Aldo Moro 2, 1-00185 Rome, Italy. 

^ Dipartimento di Fisica, Universitd di Cenova and INFN Sezione di Genova, Via 
Dodecanneso 33, I-I6I46 Genova, Italy. 

^ Dipartimento di Fisica and INFM UdR Roma 1, Universitd "La Sapienza" , 
P.zzle A. Moro 2, 1-00185 Roma, Italy. 

^ INFN - Sezione di Roma "La Sapienza". 



Abstract 

The transport of scalar quantities passively advected by velocity fields with a small- 
scale component can be modeled at meso-scale level by means of an effective drift 
and an effective diffusivity, which can be determined by means of multiple-scale 
techniques. We show that the presence of a weak large-scale flow induces interesting 
effects on the meso-scale scalar transport. In particular, it gives rise to non-isotropic 
and non-homogeneous corrections to the meso-scale drift and diffusivity. We dis- 
cuss an approximation that allows us to retain the second-order effects caused by 
the large-scale flow. This provides a rather accurate meso-scale modeling for both 
asymptotic and pre-asymptotic scalar transport properties. Numerical simulations 
in model flows are used to illustrate the importance of such large-scale effects. 

Key words: Scalar transport, transport modeling, multiple-scale analysis 
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1 Introduction 



Diverse scientific disciplines cope with systems characterized by interactions 
across a large range of physically significant length scales. This is for instance 
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the case of earthquake dynamics (see, e.g., Rcfs. [1,2]), where the important 
spatial scales range from the small scale (1 cm to 1 m), associated with friction, 
to the tectonic plate boundary scale (10^ — 10"^ km). The range of temporal 
scales goes from seconds (during dynamic ruptures) to 10^ — 10^ years (repeat 
times for earthquakes) to 10^ — 10^ years (evolution of plate boundaries). 

Many active scales arc also encountered in biological and soft matter sciences. 
For instance, to understand the physics of a cell membrane, one starts from 
Angstrom-sized atoms with motions in the femtosecond range, to go to the 
whole cell, whose diameters can be tens of micrometers and lifetimes of the 
order of days. All these scales can neither be probed by a single experimental 
technique, nor numerically resolved in a single simulation (see e.g. Ref.[3]). 

A similar situation can be found in the framework of atmospheric or ocean 
sciences. Focusing on tracer dispersion, advection by geophysical flows results 
in the rapid generation of fine-scale tracer structures. This is true even though, 
e.g. in the atmosphere, the wind field is typically dominated by large scale fea- 
tures, such as synoptic scale weather systems. Coarse-grained representations 
of atmospheric winds can indeed generate fine-scale tracer filaments, through 
the chaotic advection process (see e.g. Ref. [4]). 

A common and long-standing problem linking all these scientific disciplines is 
to build a coarse-grained description starting from the microscopic model. The 
challenge is thus to reduce the typically huge number of degrees of freedom 
by modeling the fastest and smallest scales. In this way one may construct 
a computationally tractable effective equation, which involves only the scales 
one is interest in [5,6]. 

Our aim here is to shed some light on this very general problem in the spe- 
cific framework of tracers dispersion. In this context, the scales of interest 
are generally larger than the diffusive length-scale. The problem of developing 
coarse-grained models, closed on the large scales, is also known as "parameteri- 
zation problem" [7] . For general advecting velocity fields, the parameterization 
problem is not tractable by means of systematic approaches able to deduce 
the form of the closure starting from the original, full-scale, equations. 

An important exception occurs in the presence of scale separation between 
the advecting velocity and the tracer i.e., when the former is concentrated 
at small scales and one looks at the tracer dynamics on scales much larger 
than those of the velocity. In this hmit the goal is to derive the expression of 
the asymptotic diffusion coefficient renormalized by the presence of the small 
scale velocity field. This can be accomplished exploiting asymptotic methods 
(see, e.g., [6,8,9,10,11,12,13,14] among the others). However, in many physi- 
cal circumstances one has that the velocity field may be though as a small- 
scale advecting velocity field (at scale i) superimposed to a large-scale, slowly 
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varying component (at scale L S> ^). A physical example of this situation is 
represented by the dispersion of a pollutant that, released in the atmospheric 
boundary layer (having characteristic time scales of the order of minutes), is 
observed at synoptic scales (i.e. over weeks or more). Scale separation is guar- 
anteed by the fact that meso-scale structures (in the range between weeks and 
minutes) turns out to be much less energetic than the other active scales (see 
e.g. Ref. [15]). 

In this case two levels of coarse-grained description are of interest. As in the 
previous case one would like to understand the asymptotic diffusive regime at 
scales much larger than L. On the other hand, one is also interested in deriving 
meso-scale models for the (pre-asymptotic) transport of a tracer in the range 
of scales between I and L. The latter issue has been recently addressed in 
Refs. [16,17]. In particular, explicit expressions for the effective (meso-scale) 
diffusivity and velocity have been obtained in the limit of strong large scale 
flows [17]. This investigation pointed out the important fact that the (meso- 
scale) diffusivity does depend on the large-scale advecting velocity, as well as 
(of course) on the small scales. 

Our main aim here is to quantitatively understand whether or not the opposite 
limit (i.e. a large-scale advective component which is weaker than the small- 
scale advection) leads to explicit expressions for the effective parameters which 
functionally depend on the large-scale advecting velocity. 

As we shall see, a class of eddy-diffusivity fields will emerge from a pertur- 
bative approach through which both the asymptotic and the pre-asymptotic 
transport properties can be successfully described when the effects of the large 
scales are properly taken into account. In particular, by means of numerical 
simulations of model flows we show that the meso-scale transport model is 
able not only to recover the asymptotic properties but also to predict the 
pre-asymptotic regime. The latter issue is conceptually connected with the 
predictability problem of the second kind [18,19], where the goal is to predict 
the evolution of the system through a model in which not all the degrees of 
freedom are resolved. 

For the sake of simplicity our analysis will be restricted to two-dimensional 
flows. It should be noted that working in two-dimensions is still relevant to 
many applications, e.g., to investigate the time- varying transport and mixing 
properties of isentropic flows in the atmosphere (see e.g. Ref. [20]) and in the 
ocean in connection with horizontal geostrophic eddies (see e.g. Ref. [21]). 

The material is organized as follows. In Sec. 2 the general framework for the 
multiple-scale analysis for a passive tracer in the presence of a slowly vary- 
ing advective velocity superimposed to a small-scale component is described, 
and the limit of weak large-scale advection is explicitly considered. In Sec. 3 
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eddy diffusivities are computed in two commonly considered model flows. In 
Sec. 4, the meso-scale model is compared with direct numerical simulations 
of the original transport problem both for the asymptotic and pre-asymptotic 
properties. Conclusions and perspectives are reserved to Sec. 5. The Appendix 
presents some technical material. 



2 Multiple-scale analysis 

We consider the evolution of a passive scalar field, 9{x, t), in an incompressible 
velocity field v: 

dt9{x,t) + ivd)dix,t) = DodAeix,t) , (1) 

where Dq is the molecular diffusivity. Following Refs. [16,17], we focus on 
situations where v can be thought as the superposition of a "large-scale" 
velocity field U{x, t) and a "small-scale" component u{x, t) which are assumed 
to vary on length-scales of order 0{L) and 0{i), respectively. Scale separation 
between the two fields is measured by the small parameter e = £/L. In the 
limit e — 0, multiple-scale analysis provides a description for modeling the 
dynamics of the scalar field at meso-scale, i.e. at scales larger than i and of the 
same order of L, in which the dynamical effects of the smallest scales appear 
via a renormalized (enhanced) diffusivity [16,17]. 

Following the multiple-scale approach we introduce a set of slow variables 
X = ex, T = e^t and t = et in addition to the fast variables {x,t). The 
scaling of the times T and r are suggested by physical reasons: we are searching 
for diffusive behavior on large time scales of 0(e~^) taking into account the 
effects played by the advection contribution occurring on time scales of 0{e~^). 
According to the prescription of the method, the two sets of variables are 
treated as independent. It then follows that 

di^di + eVi, dt^dt + edr + e^dr, (2) 

u^u{x,t), U^U{X,T), (3) 

where d and V denote the derivatives with respect to fast and slow space 
variables, respectively. 

The passive scalar field is expanded as a perturbative series 

e{x, t; X, T; r) = + e^^^) + e^^^') + 0{e^) . (4) 

By inserting Eqs. (2) and (4) into Eq. (1) and equating terms with equal 
powers in e one obtains a hierarchy of equations. By imposing the solvability 
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conditions on the first two orders in e one derives a Fokker-Planck equation 
for tfie "large-scale" scalar field defined as Ol = (^(°^) + e{e^^^) (see Ref. [16] 
for a detailed derivation). 

The cddy-diffusivity tensor of this "coarse-grained" Fokker-Planck equation is 
given by 

D,,{X,T)=6^,Do-{uaj) (5) 
where brackets indicates spatial and temporal averages over the fast variables 
X and t and x(a;, t; X, T) is an auxiliary field with vanishing average over the 
periodicities that obeys to the following dynamics 

dtXj + [{u +U)-d] Xj - Do d'xj = -uj . (6) 

The meso-scale transport equation for the "large-scale" scalar field 9l reads 

dtOL + {U ■ d)eL = {D^jdjOL) (7) 

where, for the sake of notation simplicity, the usual variables t and x have been 
restored. We remind that Eq. (7) is the Fokker-Plank equation corresponding 
(in the Ito convention) to the Lagrangian description 

^ = + B.,V, (8) 

where — Ui + djDij is the effective meso-scale velocity, BikBjk = Dij + Dji 
and rjiS are zero-mean Gaussian variables with {'ni{t)r]j{t')) — SijS{t — t'). 

It is worth noticing that the multiple-scale approach reduces the calculation 
of eddy diffusivities and meso-scale velocities to the solution of the auxiliary 
equation (6). In generic flows, when U is not a constant mean flow but depends 
on X and T, the equation must be solved for each value of U. This can be 
rather demanding in terms of computer resources if numerical methods are 
required to solve such an equation. Therefore, except for a few cases, in which 
analytic solutions of Eq. (6) are available (e.g. in the case of orthogonal shears 
[16], this approach does not provide a practical tool for evaluating the eddy- 
diffusivity. 

2.1 Weak large-scale flow asymptotics 

The natural way to overcome the problem of finding solutions of Eq. (6) for 
arbitrary values of the large-scale flow is to try to find explicit, even if approx- 
imate, expressions for the eddy-diffusivities. In this section we show that this 
program can be fulfilled in a perturbative limit. To be more specific, here we 
consider situations in which the intensity of the large-scale flow is weak com- 
pared to the small-scale one (the opposite hmit was considered in Ref. [17]). 
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In this case it is possible to find the solution of Eq. (6) as a perturbative 
expansion in the small parameter e — U/u: 

X(x, t; X, T) = + ex^'^ + e\^'^ + . . . . (9) 
Consistently, from Eq. (5) the eddy diffusivity can be written as 

T) = Do5ii - {uaf) - e{u,xf) - e\u,xf) + 0{e^) . (10) 

Defining the 0(1) field, U' — U /e, plugging Eq. (9) into Eq. (6) and equating 
terms with equal powers in £, one obtains the following hierarchy of equations: 

+ {u ■ d)x'-'^ ~ Dod^x^'^ = -u , (11) 

dtx^'^ + {u ■ a)x(') - Dod'x^''^ = -{u' ■ a)x(°) , (12) 



+ {u ■ a)x(") - Dod'x'--^ = -{U' ■ 9)x("-'^ . (13) 

The zeroth-order solution does not depend on the large-scale variables, and 
can be formally written as: 

x\^\x, t) = - j G{x- x', t - t')ui{x', t') dx' dt' , (14) 

where G{x,t) is the Green function associated to the linear differential oper- 
ator of the Fokker-Planck equation: 

jCfpG(x, t) = [dt + (u-d)- Dod'^] G{x, t) = 5{x, t) . (15) 

In the same way the n-th order solution can be written as 

xt'\x, t;X,T)^-Y: U'jiX, T) f G{x-x' , t-t')d^xt~'\x\ t'; X, T) dx' dt' . 

(16) 

Let us notice that all the odd terms of the expansion in Eq. (10) must vanish 
by symmetry: due to the recursive nature of Eq. (16) they correspond to 
polynomial correction of odd order in U to the eddy diffusivity, which would 
depend on the sign of U . 

Plugging Eqs. (14) and (16) into Eq. (10) we can finally write the following 
polynomial expansion in U for the eddy diffusivity: 

A,(X, T) ^D^ + Y: UiUmT\^ + 0{e') (17) 

Im 

where we denoted with Df- the eddy diffusivity originating from the small-scale 
velocity field when the large-scale one is neglected, i.e. : 

A? = D,5,i - {u,xf) . (18) 
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It is important to note that the coefficients 



= {u^G * {di{G * (9^(6- * uj))))) (19) 

(where * denotes the convolution integral) are determined only by the small- 
scale flow characteristics. We remark that this approach is relevant also in 
those cases in which direct (analytical or numerical) computation of these 
coefficients is not possible, since they can be guessed on an empirical ground. 

It is worth recalling that the expression (17) for the eddy diffusivity is a 
perturbative expansion in two different small parameters: e = £/L (which 
measures the scale separation between the two component of the flow), and 
e — U/u (which measures the ratio of their intensity) . The flrst expansion is 
the basis of the multiple-scale approach. The second allows to disentangle the 
dependence on the large-scale flow in the eddy diffusivity, since it provides 
recursive expressions for the auxiliary flelds x*^**) as polynomial expansion in 
U. 

In order to illustrate how Eq. (17) can be used to retain the effects induced 
by the large-scale ffow on the eddy diffusivity in meso-scale modeling, let us 
concentrate on the flrst non-vanishing correction, i.e. the second-order one in 
e — u/U: 5dIj\x,T) = —e'^{uixf^)- For the sake of simphcity, we consider 
the case of a two dimensional flow in the case with u is statistically isotropic. 
Without loss of generality one can choose the x'l axis parallel to U. In this 
coordinate system 

5Df],{X, T) = U\X, T){ui,G * {dv{G * {dv{G * Uj,))))) (20) 

The off-diagonal terms 6D^^\,{X, T) and 6D^2'v{-^^ vanish by isotropy. The 
correction to Df- is then diagonal in the reference frame with axis X[ || U and 



X2 -L U . Finally we can write 



A,-(Jf , T) = L>g + 5Df^{X, T) + 0{e^) , (21) 



where 

5Df^{X,T)^U\X,T)R{<t>) 

with 



^ a 

13 



R'{4>) , (22) 



1 / 



(23) 



V 



being the rotation matrix and the angle between U and Xi. Therefore the 
effects induced by the large-scale flow on the eddy diffusivity, up to the second 
order, can be obtained in terms of two parameters only. These are determined 
solely by the small-scale features: 
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a = (ii||G*(a||(G*(a||(G'*«||))))), 



(24) 
(25) 



In the general three-dimensional case the transverse correction (3 pertains to 
the plane perpendicular to the direction of the large-scale flow. 



3 Computation of the eddy diffusivity in model flows 

Let us now discuss the approximation previously obtained in two model flows. 
As representative examples of two broad classes of realistic instances we focus 
here on two large-scale flows: case (a): a steady shear, the Kolmogorov flow [22] 



where L — 2n/K is their characteristic length-scale. The two above flelds 
correspond to two typical situations: the shear flow strongly enhances the 
large-scale diffusion coefficient (in the shear direction), while the cellular flow 
(due to trapping) is characterized by a weaker enhancement of the diffusivity. 

Concerning the small-scale velocity component we consider a small scale replica 
of the cellular flow (27), i.e.: 

u — {usm{kx) cos{ky), —ucos{kx) sm{ky)) , (28) 

with characteristic length-scale given hy i = 27r/k and amplitude u. In the ab- 
sence of large-scale velocity flelds and for high Peclet numbers {Pe — u£/Dq), 
it is possible to show [24,25] that this periodic array of small vortexes give rise 
to an enhanced effective (small-scale) diffusivity ~ Do\/ Pe. A precise es- 
timation of D'^ , which is indeed in good agreement with the above expression, 
can be obtained by the numerical solution of Eq. (6), with [7 = 0. In particu- 
lar, one flnds an isotropic eddy diffusivity induced by the small-scale cellular 
flow D^j — D^Sij. The numerical technique used to solve Eq. (6) is based on 
a standard implementation of a fully de-aliased pseudospectral algorithm on 
a biperiodic square lattices of size L — 2n with 512^ collocation points. 

Let us now discuss the effect on the eddy diffusivity tensor induced by a weak 
large scale flow. 

In order to estimate the second-order correction and to assess the limits of va- 
lidity of truncating the perturbative expansion at second order in e — U/u, we 



U = {Usin{Ky),0) 
and case (6): a large-scale cellular flow [6,23,11] 



(26) 



U — {U sin(i^x) cos{Ky), —U cos{Kx) sin{Ky)) , 



(27) 
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Fig. 1. a (solid line) e P (dashed line) computed as a function of the ratio £ = U/u. 
Simulations have been performed with a small scale cellular flow (28) with u = 1 
and k = 8 (i.e. £ = 7r/4), the molecular diffusivity is Dq = 10~^ (corresponding to 

Pe = 785) and U = {U,0). 

proceeded as follows. The full auxiliary Eq. (6) is solved for various choices of 
e. As already stated, the coefficients a and f3 defined in Eq. (22) are indepen- 
dent (in the case of small-scale isotropic fiows) of the detailed spatio-temporal 
structure of the large scale velocity, which is therefore taken, in a suitable ref- 
erence frame, as a constant flow in the x-direction {U — {U, 0)). Notice that 
this is consistent with the assumption of complete scale separation, i.e. e = 0, 
in which the small-scale flow "sees" the large-scale velocity as a constant mean 
flow. 

The coefficients a and /?, are then obtained by normalizing with U"^ the eigen- 
values of Dij — Dfy As shown in Figure 1 in a wide range of e values the 
coefficient a and /5 do not depend on £, meaning that the second-order expan- 
sion holds. Of course for larger values of e higher order corrections must be 
taken into account. The results shown in Fig. 1 refer to the case of a small- 
scale cellular flow with m = 1, £ = 7r/4, = 10~^ corresponding to a large 
value of the Peclet number Pe = u£/Dq = 785. Qualitatively similar results 
have been obtained for other sets of parameters. 

However, since the small-scale cellular flow is not exactly isotropic. Therefore a 
and (3 are expected to display a weak dependence on the direction of the large- 
scale flow. To test this angular dependence, we repeated the measurements of 
a and /3 keeping flxed the ratio e = U/u (within the interval of validity of the 
second order approximation) and varying the angle between the direction 
of U and the x-axis of the small scale cellular flow. As one can see in Fig. 2a, 
the angular dependence is confirmed to be very weak. Thus, for any practical 
use, the small scale cellular flow can be considered almost isotropic. We also 
measured the angular discrepancy between the direction of U and the eigen- 
directions of SD^f as a function of (f) (see Fig. 2b). Also in this case the angular 
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Fig. 2. (a) The coefficients a e /? as a function of angle (j) that defines the direction 
of the large-scale velocity field, (b) Angular discrepancy ijj between the direction of 

the large-scale flow and the eigen-direction of the eddy diffusivity tensor. Here the 
large scale flow is chosen as U = ?7(cos((^), sin(0)), where U is fixed as i7 = O.Olii. 
The other parameters are as in Fig. 1 
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dependence is indeed rather small. Therefore, in the following we shall ignore 
it and substitute a{(f)) and (3 {(f)) with their angular averages. 

Once the correction to the small scale diffusivity tensor have been parame- 
terized through the two coefficients a and P, one can reconstruct the spatial 
structure of the effective eddy diffusivity for a generic large scale flow, within 
the second order approximation by using Eq. (22). 

In Figure 3 we show the eddy-diffusivity Dij{y) resulting from the small- 
scale cellular flow (28) superimposed to a large-scale shear (26) . We compare 
the exact solution of Eq. (6), Dij{y), with its second-order approximation for 
£ = 0.1. As one can see in Fig. 3, though we are at the border of the vahdity 
interval of the perturbation theory (see Fig. 1), the second-order approxima- 
tion recovers quite well the exact multiple-scale solution, both in the parallel 
and perpendicular direction with respect to the large-scale shear. It should be 
remarked that the improvement brought by the second-order approximation 
with respect to the zeroth-order one (the constant dotted line in Fig. 3) is 
impressive. 

It is worth noticing that the corrections induced by the large scale are not 
only non-homogeneous but also non-isotropic. In particular, the diffusion is 
enhanced in the longitudinal direction (being a > 0) while it is decreased 
in the transverse direction (/5 < 0). As we shall see in the next section this 
anisotropy is crucial to recover the pre- asymptotic features of scalar transport. 
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Fig. 3. The eddy-difFusivity Daiy) normalized with the small-scale diffusivity D^^. 
The second-order approximation (dashed line) for D^xiy) (above 1) and Dyy{y) 
(below 1) is in almost perfect agreement with the exact multiple-scale solution 
(solid line). Data refer to the case of a small-scale cellular flow (28) superimposed 
to a large-scale shear (26), with parameters values U = 0.1, L = 27r, u = 1, i = L/8, 
and Do = 10~^. The scale separation is e = 1/8 (see text for a discussion about this 
point) . 

4 Asymptotic and non-asymptotic scalar transport 

Up to now the discussion essentially focused on the formalism and the methods 
to obtain the effective meso-scale description, it is now natural to wonder 
about the practical use of this approach. As first, to test the consistency of 
the approach, we check whether the model recovers the asymptotic behavior 
of the scalar field at very large scales C ^ L and very long times T ^ T. 
Then we shall show to what extent Eq. (7) is a good model for the meso-scale 
transport, i.e. in the interval i <^ r < L. In particular, we shall show that, 
only properly taking into account the effects of the large-scale flow on the 
diffusivity tensor (i.e. at the second order Eq. (17)), the meso-scale evolution 
of the scalar field can be correctly predicted. 

4-1 Macro-dynamics of scalar transport 

As it is well known (see Ref. [6] for a review on the subject), in the asymptotics 
of large scales and long times, scalar dynamics reduce to an effective diffusion 
equation : 

dr9c = Df^,V,ec , (29) 

where 6^ is the scalar field averaged over volumes of size L. Following the 
strategy devised in Ref. [17] (see also the Appendix), we notice that the very 
large scale diffusivity tensor Df-, can be obtained in two ways. The first is 




11 




15 
10 
5 

' ' ' ' ' ' 

0.5 1 1.5 2 2.5 3 

Fig. 4. Time evolution of the asymptotic eddy-diffusivity, computed via Eqs. (A.l) 
and (A. 3), up to convergence to its constant value (non-dimensional units). The 
velocity field is given by the superimposition of a large-scale cellular flows (27) and 
a small-scale one (27), with parameters values U = 0.1, L = 27r, u = 1, £ = L/8, and 
Dq = 10~^. The time evolution of the asymptotic eddy-diffusivity (solid line) is 
well approximated by meso-scale modeling (7) in which the second-order correction 
induced by the large scale flow are retained (dashed line), while the zeroth-order 
approximation (dotted line) does not match the actual value. Time is normalized 
with the large scale advccting time scale T = L/U. 

the obvious implementation of the multiple-scale analysis directly to Eq. (1). 
On the other hand one may proceed by using two successive homogenization 
steps. In the first one wipes out the small-scale details of the velocity field, 
this is equivalent to obtaining Eq. (7). The second step consists in applying 
the multiple scale technique to the meso-scalc equation (7). In principle if 
more than two scales are present the method can be iterated (see Ref. [17] 
for a detailed discussion of this point). Both procedures lead to an effective 
diffusive equation (29) but with two (a priori) different eddy diffusivity tensors. 

The former method gives the exact value of the eddy-diffusivity tensor 
(sec Eq. (A.l)), hut requires the detailed knowledge of the velocity field at both 
large and small scales, which is a far too unrealistic situation in applications. 
The second procedure (see Eq. A. 3) is less accurate but is based on the sole 
large-scale velocity, U, and the eddy-diffusivity Dij(X,T) that parameterizes 
the small-scale flow and can be obtained with the procedure described in the 
previous section. 

In order to test whether the meso-scale model based on the perturbative ex- 
pansion (17) is able to recover the asymptotic behavior of the scalar field, we 
compare the exact value _D^'^^ obtained by a direct homogenization of Eq. (1), 
with the approximations D^''^^ and D'~'^'^. These latter are obtained by ho- 
mogenization of meso-scale model (7) where the zeroth-order approximation 
Dif^^^ — Df- and the second-order approximation D^^"^^ — Df- -|- ?>Df^ are 



12 



iype 


on 

tj Li 


11 111 


U ' 


U ' 


U ' 


s 


1 
8 


0.1 


Dn,x = 0.486 


0.520 


0.429 








Dyy = 9.53 • 10-3 


8.64 • lO^^ 


1.20 • 10^2 


c 


1 
8 


0.1 


^xx ~ ^yy ~ 3.25 • 10 


3.05 • 10-2 


3.74 • 10-2 



Table 1 

Asymptotic eddy-diffusivity resulting from the small-scale cellular flow (28) {u = 
1.0, i = 7r/4) and large-scale (S) shear (26) or (C) cellular flow (27) {U = 0.1, 
L = 27r). In all cases we used molecular diffusivity Dq = 10~^. 

used for the eddy-diffusivity (see Eqs. (22)). 

Results are summarized in Fig. 4, where we show the time evolution of the 
asymptotic eddy-diffusivity up to its convergence to a constant value in the 
case of the large-scale cellular flow (27). Notice that by using the second-order 
approximation one obtains a 6% discrepancy from the exact value with respect 
to a 15% for the zeroth-order approximation. In Table 1 one may directly 
compare the different approximations for the case of large-scale cellular and 
shear flow. 

In all the investigated cases taking into account the effect of the large scales, 
even if at the lowest nontrivial order (i.e. the second one) allows for an im- 
provement of at least a factor 2 in the relative error. Though not astonishing 
this goes in the correct direction and confirms the consistency of the approach. 
We conclude by noticing that here we are comparing "global" quantities while 
we expect a better performance of the meso-scale model for local properties 
at the intermediate scales, for which the model itself has been developed for. 

4-2 Meso-scale dynamics of scalar transport 

In many applications the asymptotic properties are less important than the 
meso-scale ones. For instance, let us consider an initially localized concentra- 
tion field (as, e.g., a pollutant released in a given region). More than being 
interested in the time-scales necessary for it to distribute uniformly in the 
whole domain, one is interested in predicting the spatial patterns and evo- 
lution of the pollutant concentration at intermediate time and scales (this is 
dramatically important if, e.g., the pollutant is a toxic substance). 

In this perspective it is interesting to see if the model (7) has a "predictive" 
character for the scalar dynamics at these intermediate scales and finite times. 

In order to test such a possibility we devised the following strategy. We con- 
sider a scalar concentration initially localized at the small scales. Then we 
follow its evolution in a square domain with periodic boundary conditions 
according to: 
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i) the exact dynamics of the scalar field ^(t) given by Eq. (1); 

ii) the model dynamics given by Eq. (7) with Dij — d\^^^ — Df-^ in the 
following this is called the ^(mo) field; 

iii) the model dynamics given by Eq. (7) with the refined approximation for 
the diffusivity tensor Dij = D^f^^ = + SD^f (see Eqs. (22)), in the 
following we denote as 9(m2) the resulting field. 

Finally the fields evolving with the three above dynamics are compared. The 
comparison is done at different levels. We looked both at the evolution of local 
and global quantities. In particular, we compare the evolution of the variance 
of the three fields, i.e. 

= {OUi^, t)) - {9(o.){x, t)f , a = T, MO, M2 , (30) 

where the average is performed over the spatial domain. The behavior of a"^^-^ (t) 
describes the mean decay of the scalar fiuctuations due to the joint effects of 
molecular dissipation and advection. Natural indicators to characterize the 
degree of spatial "similarity" between the different fields are: 

E(M2){t) = ((%)(aj,t) - e^M2){.x,t)f)/a}r^{t) . (31) 

With the normalization by crp,j(t), E(^Mo)if) and E(^Mo)if) provide a measure 
of the relative distance between true and model fields. It is worth noticing 
that, since 6'(mo) and 6'(m2) models the "true" evolution only at scales larger 
than £, in principle one should measure the errors (31) with a filtering out 
these small scales from 9{t). Here, to avoid the arbitrariness of the choice of 
filtering procedure, wc did used the unfiltered field 9(t)- This introduces an 
additional error caused by the fast-decorrelating small-scales features, which 
nevertheless becomes negligible as soon as the concentrations O^a) spreads over 
scales significantly larger than ^. 

On the numerical side, the three fields arc identically initialized as Gaussian 
distributions with width ~ 0(£) and centered in the same point of the domain. 
The robustness of the results we are going to present has been tested by 
repeating the computation with different initial locations. In the following we 
present the results for both the large-scale shear (26) and cellular flow (27). 

In Figs. 5a,b we show the time evolution of the above defined indicators. 
As one can see after a short transient the second-order meso-scale modeling 
recovers the actual decay rate, while the zeroth-order approximation does 
not. Moreover, as evidenced in the insets, the relative errors between the true 
evolving field and its two models MO and M2 is such that, while the zeroth- 
order approximation rapidly goes toward 100% error, the second-order one 
remains below 15% during the whole evolution. 
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Fig. 5. (a) Temporal decay of the fluctuations, a'^j,^ (bold line) t^^^q^ (dotted line) 
and 0"^^j2) (dashed line), in a large scale shear and a small scale cellular flow. The 
inset shows the normalized error E(^MQ-^{t) (dotted line) and E(^j^2){t) (dashed line), 
(b) The same for the case of large scale and small scale cellular flows. Parameters as 
in table 1. Time is normalized with the large scale advecting time scale T = L/U . 




Fig. 6. Snapshots of 0(^/0) (l^ft column), O^t^ (middle column), and 9(^m2) (right col- 
umn) at three different time, from top to bottom t ~ 0.5, 1.0, l.ST where T = L/U 
is the large scale advection time. The large-scale velocity field is the shear flow (26). 
At t = 9a s were initialized as a Gaussian centered in (9/16L, 7/16L) with a width 
2a ~ L Simulation parameters are summarized in table 1. 
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The improvement brought by the second-order corrections is even more strik- 
ing if one looks directly to the snapshots of the concentrations fields (Figs. 6, 7). 
As one can see the spatial patterns of the MO approximation rapidly decor- 
relate with those of the true field, that are actually well described by the 
M2 approximation. We remark crucial for the fidelity of the second-order ap- 
proximation is that the corrections to the small-scale diffusivity tensor retains 
relevant information of the anisotropic and inhomogeneous diffusive behavior 
induced by the presence of the large-scale flow, namely the enhancement of 
diffusion in the direction on U and its reduction in the transverse direction. 
This is particularly evident in the case of large-scale shear shown in Fig. 6, 
where one can see that the blob MO spreads too quickly in the y-direction 
and too slowly in the x-direction with respect to the true field. On the other 
hand, since D^^f > D^^^^ and Djf 2) < /^igfo) ^-^^^^ features are captured 
by the M2 approximation. This is even clearer in Fig. 7, where one can see 
that, differently from the M2 model, the trapping of the concentration in the 
large-scale cells is completely missed in the MO approximation. 
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5 Conclusions 



We have studied both analytically and numerically the effects of a weak ad- 
vecting velocity field at large scale on the meso-scale modehng for the trans- 
port of passive scalars. By means of multiple-scale methods we perturbatively 
computed the dependence of (pre-asymptotic) eddy diffusion tensor Dij on 
the large-scale velocity field. The corrections to D are non-homogeneous and 
non-isotropic. In particular we find an enhancement (reduction) in the longitu- 
dinal (transversal) direction of the large-scale field. The perturbative approach 
proposed here allows to develop meso-scale models retaining (at least at sec- 
ond order accuracy) these effects, which are shown to be crucial for properly 
describe the transport dynamics. 

We conclude by noticing that the new findings obtained by our approach 
seem very promising for future applications to the numerical investigation of 
large-scale transport (asymptotic and pre-asymptotic) both in the atmosphere 
and in the ocean. Indeed the present results together with those obtained in 
Ref. [17] cover two opposite limit of transport for which explicit expressions 
for the effective parameters are available. By interpolations, one may hope to 
obtain the form of these effective coefficients under general conditions. 
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A Expressions for the asymptotic diffusivity 



There are two, in principle non equivalent, ways to obtain the large-scale 
equation (29). The first one is to apply the homogenization technique to the 
exact equation (1), while the second possibility is to start from the meso- 
scale model equation (7). Following the first approach the (exact) value of the 
eddy-diffusivity tensor, D^'^^, depends on both the molecular diffusivity and 
the advection by the total velocity field v — U + u: 
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The auxiliary field x is the solution of the following equation 



dtX + {v ■ d)x - Dod'x = -V . (A.2) 

With the second way the asymptotic eddy-diffusivity tensor results from 
the combined effects of the advection given by the large-scale flow U{X,T) 
and the effective diffusion at scale i that depends on space and time through 
Dij{X,T) (see [16] for a detailed derivation), 



2 2 2 

Here the vector field % is obtained by the auxiliary equation 

dtXk + {U ■ d)xk - di{DijdjXk) = -Uk + diDik . (A.4) 

It should be noted that even though the first procedure gives the exact value 
of the eddy-diffusivity tensor D^'^^, it requires the detailed knowledge of the 
velocity field at both large and small scales. However, the expression obtained 
from Eq. (A. 3) (which, in general, does not coincide with D^''^^) has the ad- 
vantage of being based solely on the large-scale velocity, U, with the effects 
of the small-scale fiow included in the eddy-diffusivity Dij{X,T). 
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